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We consider a dynamics generated by families of maps whose invariant density depends on a 
parameter a and where a itself obeys a stochastic or periodic dynamics. For slowly varying a the 
long-term behavior of iterates is described by a suitable superposition of local invariant densities. 
We provide rigorous error estimates how good this approximation is. Our method generalizes the 
concept of superstatistics, a useful technique in nonequilibrium statistical mechanics, to maps. Our 
main example are Blaschke products, for which we provide rigorous bounds on the difference between 
Birkhoff density and the superstatistical composition. 

I. INTRODUCTION 

o' 

^^ , Dynamics often takes place in a changing environment. This means, given some control parameter a and a local 

dynamics Xn+i = faixn) generated by some mapping fa, the control parameter a itself will also slowly change in 
time. In nonequilibrium statistical mechanics, these environmental fluctuations, if taking place on a large time scale, 
are modeled by a very useful concept, so-called superstatistics. The superstatistics concept was introduced in [ij 
and has since then provided apowerful tool to describe a large variety of complex systems for which there is change 
of environmental conditions 0Ta|- The basic idea is to characterize the complex system under consideration by a 

[\ ' superposition of several statistics, one corresponding to local equilibrium statistical mechanics (on a mesoscopic level 

r) , modeled by a linear Langevin equation leading to locally Gaussian behavior) and the other one corresponding to a 
slowly varying parameter a of the system. Essential for this approach is the fact that there is sufficient time scale 
separation, i.e. the local relaxation time of the system must be much shorter than the typical time scale on which the 

'^ ■ parameter a changes. 

II In most applications in nonequilibrium statistical mechanics the varying control parameter a is the local inverse 

temperature /3 of the system, i.e. a = /3. However, in some applications beyond the immediate scope of statistical 

^~~^ ' mechanics the control parameter a can also have a different meaning. There are many interesting applications of the 
'^ , superstatistics concept to real- world problems, for example to train delay statistics [lOj, hydrodynamic turbulence |ll| 

^^ ' and cancer survival statistics [12|]. Further applications are described in [l3l - [l8j . 

pvi In this paper we want to extend the superstatistics concept to maps, which usually have invariant densities different 

f~^ from Gaussian distributions, and thus analyze this problem in a more general context. Our generalization assumes 
that the local dynamics is not anymore restricted to a linear Langevin dynamics (as in the nonequilibrium statistical 
mechanics applications), but given by an a priori arbitrary map with strong mixing properties. We will make the 
superstatistics concept mathematically rigorous by considering simple model examples of local maps where everything 
can be proved explicitly and by estimating the error terms. In our new approach described here, we allow for a priori 
arbitrary local invariant densities pa (x) and consider a dynamics given by long term iteration of a map fa with slowly 
varying a. This leads to a mixing of various invariant densities Pa{x) with different parameter a, in a way that we 
will analyze in detail in this paper. 
V^ ] If a changes on a long time scale, long as compared to the relaxation time of the local map fa, the resulting long- 

term probability distribution of iterates is closely approximated by a superposition of local invariant densities Pa{x). 
For particular examples, Blaschke products, we will indeed provide estimates of the error terms involved and prove 
how fast the Birkhoff density approaches the superstatistical composition. On the other hand, if a changes rapidly, 
then a different dynamics arises which is not properly described by a mixing of the various local invariant densities. 
Rather, in this case one has to look at fixed points of the Perron-Frobenius operator of higher iterates of composed 
maps fa with varying a. Depending on the time scale of the changes of a, there are transition scenarios between both 
cases. 

This paper is organized as follows. In section 2, we will introduce the superstatistics concept for maps. We will 
study several examples in this section to illustrate the concept and to motivate our rigorous treatment in the later 
sections. In section 3 we present a rigorous theory for Blaschke products. We will prove the existence of invariant 
measures, determine the invariant measure explicitly as a function of the parameters involved and prove our main 
result, an upper bound on the difference between Birkhoff density and the superstatistical composition. 



II. SUPERSTATISTICAL DYNAMICS OF MAPS 

Let us consider families of maps fa depending on a control parameter a. These can be a priori arbitrary maps 
in arbitrary dimensions. Later we will restrict ourselves to mixing maps and assume that an absolutely continuous 
invariant density Pa{x) exists for each value of the control parameter a. The local dynamics is 

a;n+l = fa{Xn)- (1) 

We now allow for a time dependence of a and study the long-term behavior of iterates given by 

a^n = /a„ o/a„_i O ••■/ai(a;o)- (2) 

Clearly, the problem now requires the specification of the sequence of control parameters ai, . . . ,a„ as well, at least 
in a statistical sense. One possibility is a periodic orbit of control parameters of length L. Another possibility is to 
regard the aj as random variables and to specify the properties of the corresponding stochastic process in parameter 
space. 

In general, rapidly fluctuating parameters Uj will lead to a very complicated dynamics. However, there is a significant 
simplification if the parameters aj change slowly. This is the analogue of the slowly varying temperature parameters in 
the superstatistical treatment of nonequilibrium statistical mechanics [l|, |I^ . The basic assumption of superstatistics 
is that an environmental control parameter a changes only very slowly, much slower than the local relaxation time of 
the dynamics. For maps this means that significant changes of a occur only over a large number T of iterations. In 
practice, one can model this superstatistical case as follows: One keeps ai constant for T iterations {T >> 1), then 
switches after T iterations to a new value 02, after T iterations one switches to the next values 03, and so on. 

One of the simplest examples is a period-2 orbit in the parameter space. That is, we have an alternating sequence 
ai,a2 that repeats itself, with switching between the two possible values taking place after T iterations. We are 
interested in the long-term behavior of iterates obtained for n — > 00. Possible sequences of parameters ai, 02, . . . , a^ 
of period length L could be studied equally well, with a switching to the new parameter value always taking place 
after T >> 1 iterations. Another possibility are stochastic parameter changes on the long time scale T. 

To illustrate and motivate the superstatistics concept for maps, we will now deal with three important examples of 
families of maps fa ■ 

Example 1 We take for fa the asymmetric tent map on [0, 1] given by 

with a S (0, 1). This example is somewhat trivial, because the invariant density Pa{x) is independent of a and given 
by the uniform distribution for any value of a. Hence, whatever the statistics of the varying parameter sequence 
ai, a2, . . . is, we get for the long-term distribution of iterates given by ([2]), ([3]) the uniform distribution 

p{x) = 1 (4) 

Example 2 We take for fa a map of linear Langevin type [20l . l2l| . This means fa is a 2-dimensional map given by 
a skew product of the form 

a^n+i = 9{Xn) (5) 

y.n+1 = e^""2/„+Ti/2(a;„_5) (6) 

Here g denote the average of iterates of g. It has been shown in [20| that for r — > 0, i = nr finite this deterministic 
chaotic map generates a dynamics equivalent to a linear Langevin equation [22|, provided the map g has the so- 
called (y9-mixing property [30], and regarding the initial values 0:0 S [0, 1] as a smoothly distributed random variable. 
Consequently, in this limit the variable ?/„ converges to the Ornstein-Uhlenbeck process j22i | and its stationary density 
is given by 

pM = {^^-'^''^ (7) 

The variance parameter /? of this Gaussian depends on the map g and the damping constant a. If the parameter a 
changes on a very large time scale, much larger than the local relaxation time to equilibrium, one expects for the 



long-term distribution of iterates a mixture of Gaussian distributions with different variances /3 ^. For example, a 
period 2 orbit of parameter changes yields a mixture of two Gaussians 

Generally, for more complicated parameter changes on the long time scale T, the long-term distribution of iterates 2/„ 
will be mixture of Gaussians with a suitable weight function h{l3) for r — > 0: 

p{y)^ [ d(3 h{/3)e-^^y" (9) 



This is just the usual form of superstatistics used in statistical mechanics, based on a mixture of Gaussians with 
fluctuating variance with a given weight function [l]. Thus for this example of skew products the superstatistics of 
the map fa reproduces the concept of superstatistics in nonequilibrium statitistical mechanics, based on the Langevin 
equation. In fact, the map fa can be regarded as a possible microscopic dynamics underlying the Langevin equation. 
The random forces pushing the particle left and right are in this case generated by deterministic chaotic map g 
governing the dynamics of the variable Xn- Generally it is possible to consider any i^j-mixing map here i20j . Based on 
functional limit theorems, one can prove equivalence with the Langevin equation in the limit r — )> 0. 

Example 3: Blaschke products We now want to consider further examples beyond the immediate scope of 
statistical mechanics where the invariant density of the local map is non-Gaussian but still a full analytic treatment 
is possible. Consider mappings of a complex variable z given by 

d I, 

z — 






J=l 



where |6o| — 1 and \bj\ < 1 for j — 1,2, ... ,d and d> 2. We are interested in a dynamics restricted to the unit circle 
S^ and write w € S*^ as 

u = e*^'^'^, (11) 

so that ip e [0, 1). According to Martin 27], and as established in the following sections in much more detail, the 
invariant density of / with respect to the variable u is given by 



27r \u — zop 



Here zq = /(^o) is a fixed point of /. Blaschke products usually exhibit very strong chaotic behaviour and can be 
used as the map g in equation ([5]) if an extension to a Langevin dynamics is wished for for physical reasons. The 
remarkable property of Blaschke products is that the knowledge of a fixed point zq of the map uniquely fixes the 
shape of the invariant density p", making an analytic treatment very convenient. 

Transformation of variables u ^- ip yields the invariant density p^ with respect to the variable (p as 



p^icp) = p^u) 



du 



dtp 



gi27r<p_ 2j^|2 



1-a^-b^ 



I + a^ + b'^ — 2a cos 27r(^ — 26 sin 27riy9 



(13) 
(14) 
(15) 



Here a denotes the real part and b the imaginary part of the fixed point zq, i.e. zq = a + ib. 

For zo = we get a = 6 = and hence p^{p) — 1. This is just the invariant density of a d-ary shift map on [0,1], 
noting that for bj = the Blaschke product becomes 

z — > z'' -4==> p ^- d ■ ip mod 1 (16) 

Another example would be a Blaschke product with a real fixed point zq = a. In this case the invariant density is 

^"(^) = l + a^-"2rcos2.^ - (^^) 



A particular example, taken from [23], is bi — b2 — -^ , i.e 
The fixed point condition zq — /(zo) is solved by 



m - T^T^ (18) 



zo== ^(7-V45) w0.f45898... = a (19) 

This is a unique fixed point in the unit disk D = {z\ \z\ < 1}. 

We are now in a position to explicitly do superstatistics for Blaschke products, since the invariant densities are 
known explicitly as a function of the parameters of the map. While a rigorous treatment will be worked out in detail 
in the following sections, we here just consider a simple example to illustrate the general idea. Let us consider two 
different Blaschke products, and a periodic orbit of length 2 of the parameters. For example, we may consider an 
alternating dynamics of the two maps 

/i(z) = z^ (20) 

If we iterate /2 for a long time T, then iterate /i for the same long time T, then switch back to /2, then to /i, and so 
on, the BirkhofF density will be a mixture of both invariant densities. In the ip variable we expect to get for T -^ oo 
the superstatistical result 

Pooif) = -^Piif) + ^P2{'f) (22) 

~ 2V ^ l + a2-2acos27r(^j ^ ' 

with a w 0.145898.... This will be confirmed by our rigorous treatment in the following section. If, on the other hand, 
we switch maps after each iteration step, the result will be different. In this case we need to determine the invariant 
density of the composed map /i o f2{z) (or /2 o /i(z), depending on which map is iterated first). The composed map 
is again a Blaschke product, now with d = A: 

fu{z):=fiof2{z) = j^-^^ (24) 

The fixed point condition 

^o^/i2(^o)= /'"^i^\', (25) 

yields zq ~ 0.0464774... =: c and hence the invariant density of /12 = /i o /2 is given by 

^^^(^) = l + c2-"2Ls27r^- ^'^^ 

Similarly the other composed map is also a Blaschke product with d ~ A: 

f2i{z):^l2oh{z)^ ^f\^l[, (27) 

The fixed point condition 
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^0-/21(^0)- ;/_,%2 (28) 

yields zq = C2 sa 0.0464774... 2 sa 0.215586... and hence the invariant density of /21 = /2 o /i is given by 

P^'^^^ = l + c-Vjcos2n^ - (29) 




FIG. 1: The density pi2{'^) (lower dashed hne at (p = 0) and the density P2i{'p) (upper dashed hne at (p = 0), obtained by 
switching between /i and /2 on the short time scale T — 1. The Birkhoff density ^pi2{(p) + ^p2i{(p) corresponds to the dotted 
line. In the limit T — >■ oo, the Birkhoff density will converge to the superstatistical composition pooi'fi) (solid line). 

The Birkhoff density, which describes the long-term distribution of iterates independent of the phase of the periodic 
orbit, is then given by jpi2{'^) + \p2i{^)- 

Fig. 1 shows the densities pi2 and p2i (dashed lines), the Birkhoff density \{pi2+ P21) (dotted hne) and the supersta- 
tistical composition poo (solid line) as a function of (p. Apparently, all curves are different. But the difference between 
Birkhoff density 5P12 + 5/O21 and the superstatistical composition poo will decrease if the time scale T of switching 
between parameters is increased. There is a transition scenario from the Birkhoff density to the superstatistical result 
Poo for T — )■ cx). 

In the following sections we will make our results rigorous, by proving statements for the existence and uniqueness 
of invariant measures of Blaschke products and estimating the relevant error terms of the superstatistical composition. 



III. RIGOROUS RESULTS FOR BLASCHKE PRODUCTS 



Degree d maps of the circle 

Let Aq be normalized Lebesgue measure on the unit circle S . Suppose t : S^ 



S^ is C^ and has degree d > \. 
Then its pushforward action t* : /i h- > /ior^^ on probability measures p which are absolutely continuous with respect 
to Lebesgue measure is given by the transfer operator Cr on densities. Thus t* : pAo n- Cr{p)Xo where 



i^ripWn 



E 



e»<GT-i(e»«) 



\r'ie<)\ • 



The main example which we are concerned with is when r is the restriction to S^ of a Blaschke product 

d 



B{z)^bol[j 



j=i 



-b, 



(where |6o| = 1 and \bj\ < 1 for 7 = 1,2,..., d). 



Action of Blaschke products on Poisson measures 

For any z in the unit disk, let A^ denote the harmonic measure on the unit circle whose density pz with respect to 
the normalized Lebesgue measure Aq on S^ is given by the Poisson kernel , _^ 2 ; u E S^. 



As observed in [27] , a Blaschke product pushes forward Poisson measures to Poisson measures. 

PROPOSITION: If / : D ^ D is an analytic function on the unit disk D whose extension to D is continuous and 
where the restriction r to S^ takes values in S^. Then, for all z £ D, we have 

Proof. Given a continuous function ■(/; : S^ — > R the unique extension ^ of ■0 which is continuous on D and harmonic 
in D is given by 



where Xz is the harmonic measure on S^ determined by z e D. (See |24!] Chapter 10.) 

Given / : D — )• D an analytic function on the unit disk D whose boundary value t to S^ is continuous, mapping 
S^ to S^, then for any continuous ■;/; : S"'^ — >■ R we have that "0 o / is a harmonic function with boundary value ip o t 
whence equals ip o t (by the uniqueness theorem for harmonic extensions). 

Write T* for the action on probability measures on S^ induced by t. Then for all f/' : S^ — > R continuous we have 



1 I-2TT 



SI ^TT Jo 



1 /"^^ 



^^(^^^) dO 









^Vor)(e<)dA4C) 





(-0 O T)dXz ~ 1p O t{z) 

si 

= (0o/)(z) = v;(/(z))= / ijdXf^z) ■ 

Jsi 

It follows that T* Xz = ^f(z) by uniqueness of the probability measure. D 

Consequently, if i? is a Blaschke product with a fixed point z in the unit disk then A^ is an invariant measure for 
the action (r) on S^. Clearly r expands this measure. By Walters [2^ (Theorem 18) there is only one invariant 
probability measure absolutely continuous with respect to Lebesgue and so this must be A^. The measure is mixing 
whence ergo die. 

Blaschke products which expand Lebesgue measure 



Martin [27| gives a sufficient condition for a Blaschke product 



B{z)^b,X{^^ ^ 



(where |6o| = 1 and \bj\ < 1 for j = l,2,...,d and d > 2) to expand Lebesgue measure on the unit circle. The 
sufficient condition is 

,4^1 + 1^1 



(In the special case d — 2 and 61 + &2 = this can be improved to 






> 1 



which reduces to \bj\ < 1/^/3.) 

If a Blaschke product of degree d expands some measure on the unit circle then there must be exactly d — 1 fixed 
points on S^, but since there are d + 1 fixed points on the Riemann sphere, there must be two further fixed points, 
one in each component of the complement of S'^, and both attracting (by the Schwarz lemma). 

Dynamics on the unit disk 

Any holomorphic self map of unit disk does not increase the pseudo hyperbolic distance 

\z — w\ 



d{z,w) 



|1 — wz\ 



(a metric which is invariant under Mobius transformations preserving the disk) . In the case the self map is a Blaschke 
product of degree at least two we get contraction uniform on compact subsets of the disk. 



LEMMA: For z and w in the unit disk we have 



\z-w\ ^ |2;| + |w| 



|1-- wz| - 1 + |z||w| 

Proof. We consider w fixed and z varying around a circle centre the origin and radius r < 1. Then the pseudo 
hyperbolic distance between z and w is the absolute value of 



Inverting this yields 



and the absolute value of this equals r. Hence 



1 — wz 



V + w 
1 + wv 



which is 



givmg 



Then 



Iv + wl"^ = r^\l + wv\'^ 

vv + viB + wv + WW = r (1 + wv + wv + vvww) 

(1 — r'^ww)vv + (1 — r'^){wv + wv) + ww — r^ = . 

(1 — r ww)v + (1 — r )w\ 

(1 — r'^ww) [(1 — r'^ww)vv + (1 — r'^){vw + wv)\ + (1 — r^)' 

= (1 — r ) WW — (1 — r ww){ww — r ) 



WW — 2r WW + r ww — ww + r + r w w — r ww 



WW 



r^(l — ww)"^ 



Hence 

(1 — r ww)v + (1 — r )w\ — r(l — ww) 
which the equation of a circle in v. The centre is 

(1 — r^)'w 



and the radius is 



whence the maximum of \v\ on the circle is 



1 — r^ww 



r(l — ww) 
1 — r^ww 



(l-r2)|u;|+r(l-|w;|2) 



1 — r'^\w\ 



2U„|2 



\w\ — r'^\w\ + r — r|wp 
(1 -r|M;|)(l +r\w\) 



{1 — r\w\){r + \w\) r + \w\ 



(l-r|w|)(l +r|w;|) l + r|w| 



D 



PROPOSITION: In the case i? is a Blaschke product of degree two with opposite zeros we have 



d{B{z),B{w)) 
d{z, w) 



z + w 



1 + wz 



Proof. Introduce [z,w] — i_^^ ■ Then a degree two Blaschke product can be written B{z) — 6o[-z, 6i][z, &2]- When B 
is a Blaschke product of degree d in z then the derived map 



B'^{z,w) 



[B{z),B{w)] 



turns out to be a Blaschke product of degree d — 1 in z [23| . We prove this in the case B has degree two with opposite 
zeros 6i and &2 (written ±6) obtaining the precise formula. 



[B{z),B{w)]^ 



B{z)-B{w) 
l-B{w)B{z) 



k)[z,bi][z,b2] -6o[-w, &i][w, 62] 
1 - bo[w,bi][w,h2]bo[z,hi\[z,h2 



bo 



z^-b-' w^-b^ 



1-b^z^ 1-62^,2 



1-6, 



I. / w^—b^ \ { z 
OOO (^i_bifi,^J \J 



2-b2 



60 [{Z^ - 6')(1 - 6'w2) ~ (w2 _ b^)[l „ 52^2)j /^ _ ^2^2 



(1 - &2y)2)(l _ 52^2) __ (^2 _ 52)(^2 ^52) \^1 _ ^2^2 



60 [z2 - 62 - z'^Pw^ + b^^W^ - U.2 + 62 + ^,252^2 _ 5252^2] / ^ _ ^2^2 



bo[z^ + f'b^w^-w^-b^''z^] (I 



1 + b^b'^w'^z^ — w'^z'^ — b^b^ \1 ^ b^w^ 



bQ{z — w){z + w) f 1 — b'^w^ 



(1 — 'wz){l + wz) \ 1 — b^w^ 
Hence 



,A. ...^_ [B{z),B{w)] 



B'^iz.w) 



60(2 + w) /I — b w 



>.2,T,2 



(1 + wz) \1 — 52^2 
is a degree one Blaschke product in z whose absolute value is 

z + w 



1 + wz 



D 



COROLLARY: We have a bound on the uniform contraction on compact subsets of the disk, which is given by 

d{B{z),B{w)) ^ \z\ + \w\ ^^ 
d{z, w) ~ 1 + \z\\w\ 

Write |5| for the common absolute value of 61 and 62- When |&| < l/-\/3 put 



_ l-|b|2-^(l + |6|2)(l-3|6|2) 

'^l^'l " 2|5|2 

PROPOSITION: Given a Blaschke product B with zeros ±6 satisfying \b\ < l/-\/3 and an r satisfying r|t| < r < 1 
then the closed disk Dr centre radius r is mapped inside itself by B. 

Proof. We first show that, for |&| < 1 and r < 1, a Blaschke product B with zeros ±6 maps the disk D^ inside the 
disk Dg centre radius s where 

r2 + |5|2 



l + |6|2r2 
Treating 



1 - 62^2 

as a function of z'^ (and b'^) we can apply the lemma and obtain 



|l-52z2| - l+|z2||52| 

whence the inclusion follows since s is monotone increasing in r. 
Finally the hypothesis r > r|f,| then guarantees that s < r. D 

If r satisfies the hypotheses for two Blaschke products A and B, each with opposite zeros, then an iterated function 
system on Dr with uniform contraction K — jr^ (with respect to the pseudo-hyperbolic distance on D) is created. 
This is the situation considered by Hutchinson |26i] . 
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When iterated cyclically the two maps A and B induce a Birkhoff measure 

_. m n 

m + n ^-^ ^-^ 

i=l i=\ 

Superstatistics for Blaschke products 

The map z i— > Az (from the unit disk) is continuous relative to the supremum norm on densities, as is seen in the 
following 



PROPOSITION: For z and w in the unit disk and u in the unit circle 



\-\w\^ \-\zf 



\u — w\ 



\u- z\ 



< 2 



1 + kl \ f 1 + kl 

ii-\z\rjyii-\w\y 



\z — w\ 



Proof. 



1 — |wp 1 — |zp (1 — wui)(u — z){u — z) — (1 — zz)(u — w){u — w) 



\u — wp \u — z|^ 



\u — wpju — z|^ 



(1 — ww){uu — zu ^ zu + zz) — (1 — zz){uu — wu — wu + ww) 
\u — w\'^\u — zp 

{zz — ww){\ + uu) + [{w — z)wz + {w — z)]u + [{w — z)wz + [w — z)]? 

\u — wpju — zP 



So 



2(|z|2 - |w|2) + 2ne{[{w - z)wz + (w - z)]u} 
|u — wpju — zp 



1-kP l-|z|' 



|u — wp |u — zp 



< 2- 



(|z|-|^;|)(|z| + |^>|) + |^;-z|(|^;z| + l) 
lu — wplu — zp 



-'I" "'I (l-H2)(l-|zP) 



l+lzl 



l + |w| 



\z — w\ 



Xi-\z\Y) \{i-\w\)\ 

DEFINITION: Given a Blaschke product B and a point a in the unit disk, define the error £5(0) by 

00 

£B{a) = ^£B3,a ■ 

i=i 
Similarly, given a Blaschke product A and a point /3 in the unit disk, define the error eA{P) by 

00 

i=\ 

We are interested in the case a is the attracting fixed point of A and /? is the attracting fixed point of B. 
The individual terms £30 , a &re given by the difference in Poisson densities PBUa) ~ Pfi ^-nd the individual terms 
£yi»,/3 &i'6 given by the difference in Poisson densities /Oa»(;3) ^ Pa- These Poisson differences converge to zero, in the 
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supremum norm on densities, exponentially fast as i and j tend to infinity, by the Proposition. Hence the above errors 
are well-defined. 

Given an arbitrary composition C — Ci o . . .o Ci (or word) with Ci G {A, B} define ec '■— pc ^ PCi where pc — P7 
where 7 is the attracting fixed point of the composition C . (Thus p^i = Pa '^^ Ci = A and pd = pp if Ci = B.) 



PROPOSITION: 



2^£a»oB"oA"— » — e^(/3) 



<4- 



(l + r)2 /X"+i + if"+i 
"(l-r)4 1^ 1 -K 






,A"oB"-3 



£b(q:) 



<4 



(l + r)2 /K'"+i+X"+i 
(1 - r)4 i^ 1-if 



Proof. For all i with 1 < i < ?Ti we have 



£A*oB"oA" 



PA»oB"oA" 



Pa 



and Sa^,p = PAHfi) - Pa- Hence 



£A*oS"oA" 



£a»,^ — PA'oB"oA'"-» ~ PA'(/3) 



Furthermore 



PA' 



a™-' - PA.(«| < 2-[^^|A^(S"(7)) - A\P)\ < 2 ji±^ . 2i^^+" 



■(l-r)4' ^ ^'^^ ^'"^^' - (l-r)4 

(since (i(7,/3) < 1) where 7 is the fixed point of A'^B'^. Then, for i > m, we have 

(f + r)2 (l + r)2 



(1-rr 



(l-r)4 



(since (i(/?, a) < 1). Summing over alH > 1 gives the first conclusion. A similar argument gives the second conclusion. 

D 



COROLLARY: 



y^ Pa-oB'-oA"--- + X! Pb^oA^oB^-1 - {mpa + npfi) 



exponentially fast as m and n tend to infinity. 
The quantity 



^ m n 

{ VPA'oS"oA™-i + y,PBi 



m + n 



i=l 



J=l 



is the Birkhoff density for the cycle of maps A and B. 



saW) +£5(0;) 



oA™oB"-3 J 
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